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I. INTRODUCTION 

This paper announces results for the next-to-leading- 
order (NLO) heavy-quark corrections to the axial charge 
for protons to couple to the weak neutral current 



^ q=u,c,t q=d,s,b 



(1) 



The calculation is performed by decoupling heavy quarks 
h = t,b,c sequentially, i.e. one at a time. An extension to 
simultaneous decoupling of t, b, c quarks is foreshadowed 
in our concluding remarks. 

The charge g[l' receives contributions from both light 
u, d, s and heavy c, 5, t quarks, 

2g'^^ = {Au - Ad - As) + (Ac - A6 + At) (2) 

where Aq refers to the expectation value 



{P, s\ 97^759 \p, ' 



= 2mpSfj_Aq 



for a proton of spin and mass rrip. It governs parity- 
violating effects due to Z° exchange at low energies in 
elastic lyp and Dp scattering [0, |2| or in light atoms ^ D . 
A definitive measurement of vp elastic scattering may be 
possible using the miniBooNE set-up at FNAL Q. 

Once heavy-quark corrections ^ |^ have been taken 
into account, (7^' is related (modulo the issue of S- 
function terms at x = to the flavour-singlet ax- 
ial charge, defined scale invariantly and extracted from 
polarized deep inelastic scattering: 



9^1 



= 0.2-0.35 



(3) 



The small value of this quantity has inspired vast exper- 
imental and theoretical activity to understand the spin 
structure of the proton As a result, new experiments 
are being planned to map out the spin-flavour structure 
of the proton. These include polarized proton-proton col- 
lisions at the Relativistic Heavy Ion Collider (RHIC) 
semi-inclusive polarized deep inelastic scattering, and po- 
larized ep collider studies Full NLO analyses are 



essential for a consistent interpretation of these experi- 
ments. 

Many techniques for decoupling a single heavy quark 
are available. We rely on Witten's method |]l2j, where 
the renormalization scheme is mass independent and im- 
proved Callan-Symanzik equations [ p^ can be exploited. 
In such schemes, the decoupling of heavy particles re- 
quired by the Appelquist-Carrazone theorem [p^ is not 
manifest. However, correct decoupling is ensured by ap- 
plying the matching conditions of Bernreuther and Wet- 
zel these relate coupling constant, mass and opera- 
tor normalizations before and after the decoupling of a 
heavy quark. The advantages of this approach are its 
rigour and the fact that the final results are expressed in 
terms of renormalization group (RG) invariants. These 
invariants are Witten-style running couplings 5/i, one for 
each heavy quark h = t,b, c, and axial charges for nucle- 
ons in the residual theory with three light flavours. 

We find that, when first t, then b, and finally c are 
decoupled from (0), the full NLO result is 



2g'^f^ = (Au -Ad- As) 



V(Au + Ad + As). 



(4) 



where T' is a polynomial in the running couplings ah, 



V = 



6 
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2 ^ -1 
(5) 



and {Aq)\nv denotes the scale-invariant version of Aq de- 
fined in the following way. 

Let af = (?^/47r and j3f{af) be the gluon coupling and 
beta function for MS renormalized quantum chromody- 
namics (QCD) with / fiavours and Nc = "i colours, and 
let 7/ (a/) be the gamma function for the singlet current 



("7/^75" + d-f^j^d - 



^{qkit^ibqk) 



(6) 



fe=l 



2 



A scale- invariant current (S'^s)/ is obtained when is 
multiplied by 

Up to 0{m^^) corrections, the invariant singlet charge 
(|) is given by 

gf\.^^^E,{a3){Au + Ad + As)^ 

= {Au + Ad + As).^^ (8) 

Flavour-dependent, scale- invariant axial charges Agjinv 
such as 

can then be obtained from linear combinations of (||) and 

5^^^ =Au-Ad= (Am - Arf).^^ 

g'-^^ = Am Ad - 2 As = (Am + Ad- 2As).^^^ (10) 

Here g^' = 1.267 ± 0.004 is the isotriplet axial charge 
measured in neutron beta-decay, and g^' = 0.58 ±0.03 is 
the octet charge measured independently in hyperon beta 
decay. Taking at — 0.1, ctb = 0.2 and etc = 0.35 in (||), 
we find a small heavy-quark correction factor V — —0.02, 
with LO terms dominant. 

Our results extend and make more precise the well 
known work of Collins, Wilczek and Zee |^ and Kaplan 
and Manohar |^ , where heavy-quark effective theory was 
used to estimate 5^' in leading order (LO) for sequential 
decoupling of b and b, c respectively. Our analysis is 
also influenced by a discussion of by Chetyrkin and 
Kiihn [|l6) , who considered some aspects of NLO decou- 
pling of the t quark from the neutral current and in par- 
ticular, the requirement that the result be scale invariant. 
Related work has been done on heavy-quark production 
in polarized deep inelastic scattering using the QCD par- 
ton model and in high-energy polarized jp and pp at 
NLO (18). 

The plan of this paper is as follows. Section 2 is a 
brief review of Witten's application of improved Callan- 
Symanzik equations |]l3| to the decoupling of a heavy 
quark in mass-independent renormalization schemes. In 
Section 3, we combine it with matching conditions fl^ 
to deal with next-to-leading-order (NLO) calculations in- 
volving axial- vector currents. Section 4 is then a direct 
derivation of (|^) from the formula (^ for the neutral 
current. Our concluding remarks in section 5 indicate 
the result of extending (^ to simultaneous decoupling of 
t,b,c — done not only for numerical reasons, but also to 
check that the t, b contributions cancel for mt = rrih. 

II. WITTEN'S METHOD 

In mass- independent schemes such as MS, renormal- 
ized masses behave like coupling constants. This key 
property is exploited in Witten's method. 



Let /i be the scale used to define dimensional regular- 
ization and renormalization. Then the MS scale is 

fl ^ ^V^e-''^^ , 7 = 0.5772... (11) 

We choose the same scale p, irrespective of the number 
of flavours / being considered, and so hold jl fixed as the 
heavy quarks (masses nih) decouple: 

F —t f flavours, rrih — > 00 

Also held fixed in this limit are the coupling a j and light- 
quark masses m£ j of the residual /-flavour theory, and 
all momenta p. Feynman diagrams for amplitudes 

Af ^ AF{p,fl,aF,mep,mh) (12) 

give rise to power series in m^^ modified by polynomials 
in hi{mh/p,)- We consider just the leading power Af- 

AF=AF{l + 0{l/mh)} (13) 

As rrih tends to infinity, logarithms in Af can be 
produced by any IPI (one-particle irreducible) subgraph 
which contains at least one heavy-quark propagator and 
whose divergence by power counting is at least logarith- 
mic. The effect is equivalent to shrinking all contributing 
IPI parts of each diagram to a point. This means |l^ 
that the i^-flavour amplitudes Af are the same as am- 
plitudes .4/ in the residual /-flavour theory, apart from 
TOft-dependent renormalizations of the coupling constant, 
light masses, and amplitudes: 

AF{p,p',aF,miF,mh) 

= ^ZAA'{aF,mh/p,)A'f{p,fl,af,mif) (14) 
A' 

ttf = af{aF,rnh/p.) , mif = meFD{aF,mh/ fl) (15) 

Eventually, we will have to invert (|l5|), i.e. use a/ and 
mif as dependent variables instead of ai? and rriiF, be- 
cause we hold af and m^y fixed as rrih — 00. 
For any number of fiavours / (including F), let 

/ 

be the corresponding Callan-Symanzik operator. Then 
the amplitude Af and hence its leading power Af both 
satisfy an i^-fiavour improved Callan-Symanzik equation: 

{VF+jF{aF)}AF = (17) 



In general, both jf and Z — {Zaa') are matrices. 
If we substitute ( p^ ) in ( p7| ) and change variables. 




3 



the result is an improved Callan-Symanzik equation for 
each residual amplitude, 



{Vf+jf{af)}Af = 
where the functions [fl2l [l5| 



(19) 



/?/(«/) = T^paf (20) 
Sf{af) = Vpl-ame (21) 
7/(a/) = Z-i(7F(aF) +I?f)2 (22) 

depend on a/ alone. The lack of dependence of the 
renormalization factors in (|lj) and ( |l5| ) ensures mass- 
independent rcnormalization for the residual theory. 

Although these equations hold for any f < F, their 
practical application is straightforward only when heavy 
quarks are decoupled one at a time. So we set F = f+1, 
where just one quark h is heavy. Then it is convenient 
to introduce a running coupling jl2| 

cth = ahiaFMimh/jj)) (23) 

associated with the MSf renormalized mass m,h: 



ln(mft//i) 



dx{l- 6f{x))/(5f{x) (24) 



It satisfies the constraints 

a/i(a_F,0) = Q!_F , a,i(aF,oo) = (25) 

the latter being a consequence of the asymptotic freedom 
of the F flavour theory [F ^ 16). Also, eqs. (^, ^ 
and (|2j) imply that a.h is renormalization group (RG) 
invariant: 

Vpat = (26) 
Witten's solution of (E2f) for the matrix Z is 



Z{aF,mh/^) = exp 



"/(3h,l) 



X exp - 



dx 



Z{ah, 1) 

rd 

Pfix). 



(27) 



where "ord" indicates cc-ordering of matrix integrands 
in the exponentials. Note that it is the relative scaling 
between the initial and residual theories which matters. 
For our NLO calculation, we need the formulas 



f3f{a 



2^ / 33 _ „ 
SttVT ■' 



12tt' 



■(153-19/) + 0(x*) 



7/(-) = ^/ + ^(177-2/)/ + 0(.^) 
Sfix) = -^ + 0{x') 



(28) 



where 7 f refers to the /-flavour singlet current (o ) and 
includes the three-loop term found by Larin |l9] and 
Chetyrkin and Kiihn [p^. 



III. MATCHING PROCEDURE 

Our task is to evaluate to NLO accuracy the quantities 
ail, Cif{(ih, 1) and Z{ah, 1) in (p7|), such that the answers 
depend on a/ and not ap- 

Bernreuther and Wetzel M applied the Appelquist- 



to the gluon coupling 
constant aQ° renormalized at space-like momentum Q, 



I with h 



(29) 



and compared calculations of ag" in the F = / + 1 and 
/ flavour MS theories. This reduces to a determination 
of the leading power of the one-/i-loop MSf gluon self- 
energy. The result is a matching condition 

a-^^-a-^ = CLoln— + CNLO + 0(«/,m,7i) (30) 



with a/-indepcndent LO and NLO coefficients given by 
1 



Clo — IT- , Cnlo = 
As a result, we find: 

af{ah,0) ^ ah + 0{al) = ah 



(31) 



(32) 



Bernreuther and Wetzel showed that it is possible to 
deduce all LO and NLO terms in ( |30| ) from ( pl| ) and /?/ 
and ^/ in (p^. We have done the calculation explicitly: 



NT,n I 



— In h c f In 

37r p ^ 



+ df In 
142 - 19/ 



1 



1 + ^ln^ 
rrih 



In- 



Stt V 2 / 

57 + 16/ 



27r(31-2/) ' 27r(33-2/)(31-2/) 
From (E4), we have also found oih in NLO, 



(33) 



NLO 



'-3^(1-/)'" 



153 - 19/ 
27r(33-2/) 



JTT V 2 / /i 



where nih is Witten's RG invariant mass: 

rah 

fhh = rah exp I dx 5f{x) / (3f{x) 



(34) 



(35) 



If desired, hi{mh/fJ.) can be eliminated by substituting 
In ^ — In — 

LO 



31-2/ 



■In 



3tt\2 ■'J 



(36) 



Therefore the asymptotic formula for ah as mh — s- cx) is 

5, ^ 3vr/{(f -/)ln!|i-|-fc,lnln!^-fO(l)} 

^ 3(153-19/) _ 6(33 - 2/) ^^^^ 



2(33 - 2/) 



31-2/ 
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To find the matrix Z{ah, 1) in NLO, we need a match- 
ing condition for the MS amphtude F^s for /i 7^ 75/1 to 
couple to a hght quark £. We have calculated the leading 
power due to the two-loop diagram <^)X[^ : 



and expanding in at, keeping all quadratic terms: 

6167 at' 



3312 TV 



O 



(5?)} 



f/iS — 



Consequently, there is a NLO term af^/Sir'^ in Z{ah, 1) 
for h'jfj^j^h to produce ^7^75^ as mh — > 00. 



IV. HEAVY QUARKS DECOUPLED FROM J^s 

Let us adopt the shorthand notation qf for MS currents 
(97^75?) / /-flavour theory, e.g. the neutral current 

J^g' and the scale-invariant singlet current (S/is)/: 



i(w d-^s-^c- 46)5 + 0(l/mt) (47) 
5 

Next we decouple the b quark. Here, it is natural to 
define five-flavour quantities and rhb^ analogous to 
the six-flavour running coupling at and mass frit for the 
top quark: 



In 



dx 



In 



■mis 
rriis 



(48) 



= -(t-b + c- s + u 
2 ^ 



(39) 

Sf^Ef{af)iu + d + s + ...)f (40) 

We begin by decoupling the t quark. Because of 

{c - s + u - d)e = (c - s + u - d)5 + 0{l/mt) (41) 

we see that (|27| ) is non-trivial only for 

{t - 6)6 = Zq^5{u + d + s + c + 6)5 

+ i(w + d + s + c-46)5 + 0(l/mt) (42) 

5 

Since [t — 6)g is scale invariant, we have 7^? = in (p7[): 



Eqs. ( pO| ) and (21) imply that 0:^5 and mf,5 are both 
RG/=5 and RG/=6 invariant 

2?5afc5 = Peafcs , V^rhb^ ^ ^ Vefhb^ (49) 

and hence physically significant in the original six-flavour 
theory. So we write ctb and rfib for and m^s. 

Consider decoupling the b quark from (^7|). The NLO 
matching condition (pq) becomes 



^ (in ^ + (u + + s + c)4 + 0(a^, mfc5 i) 

(50) 



(•at 

26^5(^6, "^t/M) = 26^5(5(,l)exp- / 



75(3:) 

/35(x) 

(43) 



so the non-singlet current in (47) can be written 
{u + d + s + c~ 46)5 

= {l~^}E^\ab)S4 + 0{almb^^) (51) 

For the singlet current ^5 in (^^. we find 



The operator matching condition corresponds to S5 = E5{ab)^l + ^}£^4 ^^84 + ©(ag, m^g ^) (52) 



te = ^(ln!I^ + l){u + d+s + c + b)5+ Oialmt') 



mt 



and so we conclude: 



Z6^5(5t,l) = -^ + ^5? + 0(5?) 



(44) 



(45) 



taking into account the definitions (|^) and (^0|). Then 
we expand ( pT| ) and ( p2| ) in Sf,, keeping quadratic terms: 



+ 0{alb,m-l) 



125663 . 
828007r 



6167 ~ \„ 
3312;^"*r^ 

(53) 



Eq. ( ^ ) is to be expanded about at ^ with a5 held 
fixed. In that limit, the exponential tends to the con- 
stant factor i?5(Q;5) of (^. This factor combines with the 
singlet current in ( ^ ) to form the scale-invariant opera- 
tor 5*5, as required by RG/=5 invariance. The full NLO 
result is then obtained by writing 



The same technique can be applied to decouple the c 
quark from S4 in (p3) and {c — s + u — d)4 (the result of 
decoupling b from (|4lj)). That yields the final results (^) 
and (|) given in the Introduction. 



NLO 



2^6-^5 (at, 1) '^^pj^y 

+ -{u + d + s + c-ib)5 
5 



75(3:) 
/35(x) 



REMARKS 



Our results depend on two key features: 



(46) 



1. Like previous workers in this area, we decouple 
heavy quarks sequentially, i.e. one at a time. 



5 



2. Our running couplings at, on, and 5c, which corre- 
spond to Witten's prescription are all renor- 
malization group invariant. 

The restriction to sequential decoupling is numerically 
reasonable for the t quark, but dubious for the h and 
c quarks, because it amounts to an assumption that 
ln(TOc//2) is negligible compared with ln(mf,//i). This 
inhibits detailed comparison of NLO results with data, 
which ought to be carried out with NLO accuracy |^ . 

There is also a theoretical issue here: one would like 
to check that, in the limit mt = rrib, the t and h con- 
tributions cancel. However, that is outside the region of 
validity ln(mt//i) ^ hi{rrn,/p) for sequential decoupling. 

For these reasons, we have extended our analysis to 
the case of simultaneous decoupling, where the mass log- 
arithms are allowed to grow large together: 

ln(mc//l) ^ \n{mb/p,) ~ \n{mt/fl) large 

This requires a considerable theoretical development 
of matching conditions and the renormalization group, 
which we will present separately. It involves the construc- 
tion of running couplings at. Off,, etc with the following 
properties: 

1. They are renormalization group invariant. 

2. They are defined for mt ^ mb ^ mc, and can have 
a non-trivial dependence on more than one heavy- 
quark mass. 

3. In the special case of sequential decoupling, they 
agree with at, ab and Sc to NLO. 



4. For the case of equal masses, they coincide, e.g. 

at — ab for mt = mb (54) 



Then we find that the result for the simultaneous de- 
coupling of the t, b, c quarks from the neutral current is 
of the same form (Q) as the sequential answer, but with 
the sequential running couplings in (|^) replaced by our 
simultaneous couplings at, ab, and ac'. 



125663 
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6167 
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at — 



22 1 
75^"=/ 



(55) 



Notice the factorization of the terms depending on at 
and ab- Given (|54|), the factor ab — at ensures that all 
contributions from b and t quarks cancel (as they should) 
for mt = mb- 
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